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ABSTRACT 

We present a new mid-infrared (12.5 fim) map of the BN/KL high-mass star-forming complex in 
Orion using the LWS instrument at Keck I. Despite poor weather we achieved nearly diffraction-limited 
images (FWHM = 0.38") over a roughly 25"x25" region centered on IRc2 down to a flux limit of ~ 250 
mjy. Many of the known infrared (IR) sources in the region break up into smaller sub-components. We 
have also detected 6 new mid-IR sources. Nearly all of the sources are resolved in our mosaic. The 
near-IR source "n" is slightly elongated in the mid-IR along a NW-SE axis and perfectly bisects the 
double-peaked radio source "L". Source n has been identified as a candidate for powering the large IR 
luminosity of the BN/KL region (L — 10 5 Lq). We postulate that the 12 /im emission arises in a 
circumstellar disk surrounding source n. The morphology of the mid-IR emission and the Orion "hot 
core" (as seen in NH3 emission), along with the location of water and OH masers, is very suggestive of a 
bipolar cavity centered on source n and aligned with the rotation axis of the hypothetical circumstellar 
disk. IRc2, once thought to be the dominant energy source for the BN/KL region, clearly breaks into 4 
sub-sources in our mosaic, as seen previously at 3.8 - 5.0 /im. The anti-correlation of mid-IR emission 
and NH3 emission from the nearby hot core indicates that the IRc2 sources are roughly coincident (or 
behind) the dense hot core. The nature of IRc2 is not clear: neither self-luminous sources (embedded 
protostars) nor external heating by source I can be definitively ruled out. We also report the discovery 
of a new arc-like feature SW of the BN object, and some curious morphology surrounding near-IR source 
"t". 

Subject headings: Stars: formation — stars: early-type — ISM: outflows — ISM: individual (BN/KL, 
OMC 1) — infrared: ISM 



1. Introduction 

Because of its proximity (450 pc, Genzel & 
Stutzki 1989), the OMC 1 region, which in- 
cludes the Becklin-Neugebauer object and the 
Kleinmann-Low nebula (BN/KL) and the Orion 
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"hot core" , is among the most studied regions of 
high-mass star formation in the sky. An outstand- 
ing uncertainty about this region is the source of 
the ~ 10 Lq of infrared (IR) luminosity from 
BN/KL. Though many individual peaks of ther- 
mal infrared (7 to 24 fim) emission have been 
identified (Lonsdale et al. 1982), many may not 
be self-luminous. The VLA studies by Menten & 
Rcid (1995) indicate that, in addition to the BN 
object, the region contains at least two additional 
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ultra-compact radio sources, I and L, separated 
by about 3" from each other (1,500 AU in projec- 
tion). Radio source L is coincident with IR source 
"n" (Lonsdale et al. 1982). Recent mid-IR studies 
of the region suggest that the hot stars powering 
these compact H II regions may be sufficient to 
drive the enormous IR flux from BN/KL (Gezari 
et al. 1998). 

The Orion "hot core" is a warm, dense cloud 
(T > 220 K, n(H 2 ) = 5 x 10 7 cm" 3 ) in the SE 
region of BN/KL, first observed in ammonia emis- 
sion (Morris et al. 1980). An unusually high num- 
ber of hydrogen-saturated molecules are observed 
toward the hot core, perhaps due to the evap- 
oration of icy grain mantles (Blake et al. 1987; 
Brown et al. 1988). The velocity centroid of the 
"shell masers" masers surrounding radio source I is 
nearly identical to the LSR velocity of the hot core 
( v isr = 5km s _1 ), which suggests that the two are 
physically related. The absence of any IR emission 
and the morphology of the ammonia emission in 
the vicinity also suggests that radio source I is em- 
bedded in the hot core (Gezari et al. 1998; Wilson 
et al. 2000). Hence source I is probably responsi- 
ble for most of the heating and chemistry in the 
hot core. 

Sources I and n are both associated with strong 
OH, H2O, and SiO mascr emission (Johnston 
et al. 1989; Genzel et al. 1981). An expanding 
arcminute-scale complex of high velocity (v = 30 
to 100 km s _1 ) H 2 masers surrounds the en- 
tire OMC 1 region, with a center of expansion 
that coincides, within several arcseconds, with 
source n (Genzel et al. 1981; Menten & Reid 
1995). In addition to this high- velocity flow, 
a low-velocity arcminute-scale cluster of much 
brighter maser spots is associated with this ex- 
pansion (v 'expansion = 18±2 km s _1 ). As shown 
by Gaume et al. (1998), the so-called 22-GHz H 2 
"shell" masers, which were mostly resolved-out in 
the VLBI observations of Genzel et al. (1981), 
are concentrated in a 2" by 0.5" strip centered 
on source I and are oriented roughly orthogonal 
to the bipolar CO outflow emerging from OMC 
1 (Chernin & Wright 1996). Bright SiO maser 
emission lies within 0.1" of source I, and consists 
of four clusters of maser spots having velocities 
very similar to those of the H 2 shell masers. 
Greenhill et al. (1998) and Doelcman et al. (1999) 
found that these masers are concentrated into four 



linear chains located north, east, south, and west 
of source I. Contrary to prior conclusions, a new 
analysis by Greenhill et al. (2003) suggests that 
these masers trace a disk seen roughly edge-on 
with a projected rotation axis at PA ~ 45° and 
a bipolar wind impacting the surrounding dense 
medium. 

Radio CO line emission indicates that there is 
also a fast (30 to 100 km s _1 ), weakly collimated 
bipolar outflow emerging from the OMC 1 region 
with a blue-shifted lobe toward the northwest. 
The OMC 1 outflow has a mass of about 10 M Q 
and a kinetic energy of about 4 x 10 47 ergs (Kwan 
& Scovillc 1976). There are also "fingers" of H 2 
emission, presumably from shocked gas, pointing 
away from the BN/KL region to the NW and SE. 
The OMC 1 outflow and the H 2 fingers (Allen & 
Burton 1993; Stolovy et al. 1998; Schultz et al. 
1999; Kaifu et al. 2000) appear to indicate that 
a powerful ejection occurred in OMC 1 relatively 
recently. The H 2 finger system consists of over a 
hundred individual bow shocks which delineate a 
relatively wide-angle bipolar outflow toward PA w 
315° with an opening angle of more than 1 radian 
in each lobe. The proper motion vector field of the 
H 2 line-emitting knots, as observed using Hubble 
Space Telescope (HST), has been interpreted as 
an explosion that occurred about 1010±140 years 
ago (Lee & Burton 2000; Doi et al. 2002). This 
is consistent with the OMC 1 outflow observed in 
CO. 

In this paper we present the highest spatial res- 
olution map at 12.5 fim of the BN/KL region to 
date. 1 In the next section we describe our observa- 
tions using the LWS instrument on Keck I and the 
data analysis applied to produce the mosaic. In 
Section 3 we present our results, including 6 new 
mid-IR sources, and new morphologies for IRc2 
and IR source n. In Section 4 we discuss the in- 
terpretation of some of our results, especially for 
sources n, I, and IRc2. We also try to rational- 
ize the various outflows in the region and some 
curious features surrounding near-IR source "t". 
Summary and conclusions are given in Section 5. 

Note added during review: Another dataset obtained with 
LWS on the Keck I telescope has just been published by 
Greenhill et al. (2004) and covers some of the same fields 
in BN/KL as presented here. 
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2. Observations &: Data Analysis 

Observations of the BN/KL region were made 
using the Keck Observatory facility mid-IR cam- 
era LWS on UT 16 November 2002. LWS is a mid- 
IR imaging and spectroscopy instrument mounted 
on the forward Cassegrain focus of Keck I, em- 
ploying a Boeing 128 x 128 As:Si BIB array with 
a 10.2 x 10.2" field of view (Jones & Puetter 
1993). Weather conditions were poor on UT 16 
November 2002: the thermal background and at- 
mospheric transmission varied by ~ 50% through- 
out the first half of the night. Approximately 30 
individual observations (frames) were made of the 
BN/KL region at slightly overlapping positions us- 
ing the 12.5/im filter (12 - 13 (im bandpass with 
> 80% transmission). The chopping secondary 
mirror was driven at 2 Hz with a 30" E-W throw. 
Each frame was observed using the standard mid- 
IR chop-nod technique with two chopping posi- 
tions "+" (on-source) and "-" (off-source); after 
chopping with the source in chop-beam " +" , the 
telescope was nodded along the chop axis so that 
the object would sit in chop-beam and chop- 
ping would continue. Each frame was observed for 
one complete chop-nod cycle, yielding a total on- 
source integration time of 27.6 seconds per mosaic 
frame. The standard stars Peg, And, and 
Gem were also observed throughout the night for 
PSF determination and flux calibration. Our data 
are nearly diffraction limited at 12.5 /mi with PSF 
FWHM = 0.38"and a Strehl ratio of 35%. 

The LWS instrument saves all coadded data in 
individual chop and nod beams for each observa- 
tion frame. The data were reduced by first finding 
and eliminating bad pixels in the raw data and 
then differencing the chop beams for each nod po- 
sition. There were very few contaminating sources 
in our off-source chop beams: in these cases we 
merely threw out the offending chop pair and used 
the chop pair from the other nod position. Typi- 
cally there was some residual sky signal left over 
after differencing. In order to create a smooth mo- 
saic, we opted to eliminate this residual by calcu- 
lating a sky mode and subtracting the value from 
the differenced images. While this effectively re- 
duces the background to near zero, it also elim- 
inates any diffuse nebular emission from the re- 
gion. A gain map for the detector was created 
using an average of the off-source chop positions 



for each chop-nod cycle, correcting for the dark 
current which we measured at the start of the 
night. This gainmap was then divided into the 
differenced chop-beam data. We found that, after 
nodding, the telescope did not always return the 
target to the exact same position on the detector, 
so we registered the nod position data using the 
cross-correlation of the two images and then coad- 
ded to create a final image for each telescope point- 
ing. Reduced images for each pointing were then 
combined into a mosaic by registering on sources 
common to successive frames. The final mosaic is 
shown in Figure 1. Because of our frame overlap- 
ping strategy, we inadvertently missed one small 
region of BN/KL inbetween IRc2 and BN (shown 
as a dashed box in Figure 1). 

A linear astromctrical solution for the com- 
pleted mosaic was determined using the coordi- 
nates for BN (Mcnten & Reid 1995) and IRc7 
(Gezari et al. 1998). The plate scale was deter- 
mined to be 0.0805 " pixel -1 . Astromctrical errors 
were determined from the positional errors for BN 
and IRc7 (Gezari et al. 1998), the centroiding ac- 
curacy (< 0.5 pixel), and the offset from IRc7 or 
BN (whichever is closer). Positional errors range 
from 0.1" near IRc7 and BN to 0.3" at the edge of 
the mosaic. No attempt has been made to measure 
or correct for optical distortion, which is assumed 
to be small over the (very narrow) field of view. 

Due to highly variable atmospheric transmis- 
sion, flux calibration using the observed standards 
was not reliable. Instead we calibrated the entire 
mosaic using known fluxes for BN and IRc2 from 
(Gezari et al. 1998). Since the instrumental re- 
sponse varied from frame to frame in the mosaic, 
the flux calibration for any frame that does not 
contain BN or IRc2 can deviate significantly from 
those frames that do. The regions around IRc2, 
IRc7, and BN were observed multiple times, al- 
lowing us to determine a maximum variation (2a) 
in the flux solution from frame to frame of ~ 25%. 
The uncertainty in the known fluxes of BN and 
IRc2 at 12 fim is roughly 10% (Gezari et al. 1998), 
so the absolute flux calibration error (« 2a) for 
our mosaic ranges between 10 and ~ 27%. The 3a 
flux limit for our data is approximately 250 mJy 
for a compact source measured in a 15-pixel-radius 
aperture. 
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3. Results 

The 12.5 /im mosaic of BN/KL presented here 
is the highest spatial resolution mid-IR map of 
the region to date. Even at this resolution (PSF 
FWHM = 0.38") all the detected sources are re- 
solved - there are no point sources down to our 
flux limit. Many known IR sources do break up 
into multiple compact sources, but even these sub- 
components are resolved. We have also detected 6 
new mid-IR sources, one near BN which we have 
called the "southwest arc" , and 5 more which we 
have named IRcl8 through 22, following Gezari 
et al. (1998). Positions for each of the sources rel- 
ative to BN are given in Table 1 . The uncertainty 
in the relative position is dominated by the off- 
set from IRc7 or BN, as discussed in Section 2, 
and ranges from 0.1" to 0.3". We have also mea- 
sured the flux for a handful of the most compact 
sources, using a large circular aperture and an 
outer annulus to determine the background flux 
density and its variance (Table 1). The flux uncer- 
tainty is dominated by variations caused by poor 
weather conditions (Section 2) but is not larger 
than w 27%. 

The mid-IR emission from the BN/KL region is 
generally thought to fall into two categories: Rel- 
atively dense clouds heated by external UV ra- 
diation from local high-mass stars embedded in 
the region, and dense regions heated by embedded 
YSOs. Most of the mid-IR sources (except for BN 
and n) are not seen at wavelengths less than a few 
microns owing to the significant (but patchy) fore- 
ground extinction associated with OMC 1. IRc3, 

4, 5, and 6 have relatively cool mid-IR colors and 
polarizations consistent with external illumination 
(Gezari et al. 1998; Dougados et al. 1993; Minchin 
et al. 1991), most likely from sources I and n. In 
our data, IRc3 and 6 break up into smaller frag- 
ments, but there is no evidence of point sources 
which might be interpreted as embedded young 
stars. BN, source n, and IRc2, 7, 11, and 12 
all have compact quasi-elliptical morphologies sug- 
gestive of dense clumps with embedded YSOs. In 
addition, these sources have much warmer mid- 
IR colors indicative of embedded objects (though 
the clumps would also be heated by external UV 
radiation) (Gezari et al. 1998). 

The structure of IRc2 at 12 /j,m is similar to that 
seen at 4 - 5 /xm, but not identical (Figure 2). Four 



unresolved sources (A - D) have been detected 
at L'- with marginally resolved counterparts at 
4.0 and 5.0 /im(Dougados et al. 1993). The mor- 
phology of IRc2 is clearly wavelength dependent 
from 3.8 to 12.5 ^m: the eastern peak of IRc2 at 
12.5 fim lies between A and B; and the peak of 
12.5 fim emission near source C is offset from that 
source by about 0.5" to the SE. The faint emis- 
sion just north of source D, first identified as a 
"jet-like" feature by Dougados et al. (1993), shows 
up as a distinct compact source at 12.5 /im which 
we have named "E" (Figure 3). Note that none 
of these IRc2 sources have been detected at 2 /xm 
(e.g., Lonsdale et al. 1982: Schultz et al. 1999). 
Though very close (within 0.5"), source A is not 
coincident with source I. The nature of IRc2 and 
radio source I will be discussed in Section 4. 

IR source n coincides with the double radio con- 
tinuum source "L" (N and S), which has been iden- 
tified as a possible source of significant ionizing 
radiation in the BN/KL region (Menten & Reid 
1995). Source n appears as an unresolved point 
source throughout the ncar-IR (Lonsdale et al. 
1982; Dougados et al. 1993). At 12 iim, how- 
ever, source n appears elongated along a NW-SE 
position angle (Figure 3), perfectly bisecting the 
northern and southern components of radio source 
L (Figure 1) - suggestive of a circumstellar disk 
straddled by ionized bipolar cavities (or lobes, see 
Section 4). The uncertainty in the 12 ttm posi- 
tion of source n (0.1") is smaller than the ~ 0.2" 
separation between n and both components of L. 

Though limited to one bandpass, the increased 
sensitivity of our dataset over previous efforts has 
revealed some new mid-IR sources in the BN/KL 
region: IRcl8 - 22 (following Gezari et al. (1998)) 
and the "BN SW Arc". The concentric relation- 
ship to BN suggests that the SW Arc is a nearby 
shell of heated gas and dust. Due to the inhomo- 
geneous nature of the BN/KL region and the lack 
of data at other wavelengths, we can only specu- 
late on the nature of the SW arc: it may be ex- 
ternally heated by UV radiation or is possibly a 
portion of a shocked, compressed shell created by 
an outflow or jet from BN. Source IRcl8 is com- 
pact but still resolved, and may have a counterpart 
at 4 itm(Dougados et al. 1993). IRcl9 is an irreg- 
ularly shaped source just southwest of IR source 
"q" (Lonsdale ct al. 1982); it is not clear whether 
the two are related. Source IRc20 appears as a 
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"bridge" between IRc7 and 3. This feature is also 
seen clearly at 4 /im (Dougados et al. 1993). IRc21 
is just north of IRc2 and may be associated with 
a bright knot of H2 emission immediately to the 
west (see Figure 4). IRc22 appears to be associ- 
ated with IR source "t" (Lonsdale et al. 1982) and 
will be discussed further below. 

4. Discussion 
4.1. IR Source n 

There is a significant anticorrelation between 
the 12.5 /im emission in our map and NH 3 emis- 
sion from the "hot core" (Figure 5), suggesting 
that most of the mid-IR sources are at the same 
distance or behind the hot core. The very dark 
regions in our mid-IR mosaic - the roughly trian- 
gular region NW of n and the larger region SE of n 
- coincide very nicely with the NH3 emission peaks 
(the correlation is most dramatic in the 6.79 km 
s _1 channel of NH 3 ). Source n sits precisely at the 
waist of this emission. The mid-IR and NH3 con- 
tours immediately surrounding source n are highly 
suggestive of bipolar cavities, with opening angles 
of ~ 90° , aligned nearly perpendicular to the long 
axis of source n. In addition, most of the OH and 
H2O masers in the entire BN/KL region can be 
found within an extension of the conical regions 
defined by these cavities (Figure 5) suggesting that 
these regions have been cleared by an outflow(s) 
from source n and that this outflow(s) is responsi- 
ble for the masers. 

The bipolar cavity defined by the mid-IR, NH 3 , 
and maser emission, along with the elongated 
structure of n at 12 /im, strongly suggests that 
n has a circumstellar disk with a rotation axis ori- 
ented NE SW. The unresolved near-IR flux from 
n would be from the stellar photosphere and/or 
hot inner portion of the disk. The 12 /im emis- 
sion presumably arises in the outer portions of 
the disk (at hundreds of AU) , where temperatures 
are a few hundred Kelvin, and hence could be re- 
solved in our data. The FWHM at 12 /xm along 
source n's long axis (PA = 131°) is 0.57". Source 
n is unresolved along the short axis of mid-IR 
emission. If we assume a distance of 450 pc to 
OMC 1 , source n must be 255 AU in diameter and 
< 167 AU in projected height, leading to a disk 
inclination of > 50°. The disk cannot be edge- 
on, however, or source n would not be detected 



at near-IR wavelengths owing to extinction of the 
stellar photosphere and inner disk by the outer, 
flared regions. The polarization angle and linear 
fraction for source n at 3.8 /xm (Minchin et al. 
1991; Dougados et al. 1993) is consistent with this 
disk interpretation. The elongated radio source 
L, bisected perfectly by the mid-IR disk, may be 
emission from an ultra-compact bi-polar H II re- 
gion above and below the disk or perhaps a jet and 
counter-jet. The rotation axis of the disk would be 
along PA w 41°, in near perfect alignment with 
source L (N and S) and the bipolar cavities to the 
NE and SW. Because of its significant radio con- 
tinuum emission, source n has already been iden- 
tified as a possible contributor to the energetics of 
the BN/KL region (Mcntcn & Rcid 1995). Our 
data, coupled with the NH 3 observations of the 
hot core and the maser spots observed throughout 
BN/KL, suggest a disk and outflow morphology, 
as would be expected for a young, high-mass star. 

4.2. Source I and IRc2 

The strong radio continuum from source I, in- 
dicative of an ultra-compact H II region, suggests 
that it also is a significant energy source for the 
BN/KL region (Menten & Reid 1995). Source I 
has not been detected in the IR, however, indicat- 
ing that if it is indeed a hot, young star, it must 
be obscured by over 60 visual magnitudes of dust, 
placing it within or on the far side of the hot-core 
(Gezari et al. 1998). Observations of the 86 GHz 
SiO masers led to an interpretaion, at first, of an 
outflow from source I aligned with PA ~ —45° 
(Greenhill et al. 1998). New VLBA observations, 
however, suggest that the SiO masers arise in a 
disk seen nearly edge-on with a rotation axis at 
PA - 45° (Greenhill et al. 2003). In this picture, 
the water masers centered on source I (the "shell" 
masers) are due to a wide angle bipolar outflow 
from I impacting the surrounding dense natal ma- 
terial (Greenhill et al. 2003). Indeed, the shell 
masers trace very closely the NH 3 contours sur- 
rounding source I and the mid-IR contours in our 
mosaic (Figure 6), suggesting that they arise in 
the walls of a cavity excavated out of the hot core 
gas by bipolar outflows from source I (Wilson et al. 
2000). 

The nature of the IRc2 sources (A - E) is not 
entirely clear. Because of their proximity to source 
I on the sky and the fact that the mid-IR emis- 
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sion comes from a location which is at or behind 
the hot core (in which source I is most likely em- 
bedded), it seems that the IRc2 sources are asso- 
ciated with source I. The silicate extinction fac- 
tor (9.8 /im line-to-continuum ratio) for IRc2 is 
~ 50 - higher than any other source in BN/KL 
(Gezari et al. 1998). Due to the correlation with 
source structure, the extincting grains are prob- 
ably local to IRc2 itself. The dereddened flux 
densities in the mid-IR for IRc2 indicate a color 
temperature of w 240 K and a total luminosity of 
L = 1000 ± 500L Q (Gezari et al. 1998). At least 
part of the strong silicate absorption could be due 
to the hot core, which clearly overlaps some of the 
mid-IR emission (Figure 6). 

The IRc2 sources also show significant polariza- 
tion in the near-IR (Minchin et al. 1991; Dougados 
et al. 1993). The high-resolution L' polarimetry 
map by Dougados ct al. (1993) clearly shows lin- 
ear polarizations of 10 15% for the IRc2 sources 

with polarization angles roughly perpendicular to 
the direction of source I. This naturally leads to 
the interpretation that some of near-IR emission 
from IRc2 is due to scattering of radiation from 
source I. 

The wavelength dependent morphology of IRc2 
(Fig. 2) displays two important characteristics: 

1. the resolved structure of IRc2 at 12.5 ^m 
becomes more unresolved at shorter wave- 
lengths, culminating in the point-like nature 
of sources A - D at L'; and 

2. the 3.8 - 5.0 /im sources A - D generally ap- 
pear at the periphery of the 12.5 /xm emis- 
sion. 

How to interpret these characteristics, however, 
is not at all clear: the first characteristic sug- 
gests that these sources are embedded protostars 
(as indicated by the mid-IR color of IRc2); while 
the second is indicative of an externally heated 
clump of gas and dust (as suggested by the near- 
IR polarimetry results). In addition, IRc2 ap- 
pears to be partially extincted by a tongue of 
dense foreground material which is revealed by 
a small bay in the extended 12 ^m emission at 
(Aa, AS) = (5.25,-7.80), and coincides with a 
small protrusion in the NH3 emission (Figure 6), 
so patchy extinction plays a strong role in the ob- 
served morphology as well. 



If sources A - E were indeed embedded pro- 
tostars, then the total IR luminosity for IRc2, 
L = 1000 ± 500 L Q indicates that their masses 
would likely be in the 3-8 M Q range. This 
would imply that, in the 500-AU-diameter region 
where the IRc2 cluster sources are located, the vol- 
ume averaged density would be larger than 10 10 
cm~ 3 This would represent a compression by at 
least a factor of ~50 in all three dimensions rela- 
tive to the density of a typical star-forming cloud 
core, 10 5 cm~ 3 and the original cloud core would 
have to have been larger than the whole BN/KL 
region. Therefore, while it is probably not impos- 
sible to get a group of massive stars forming so 
close to each other (the challenge is avoiding the 
centrifugal barrier during the collapse), it seems 
far more likely that sources A - E are patches of 
reflecting dust rather than embedded protostars. 

In addition, the near-IR polarization and 
wavelength-dependent morphology of IRc2 from 
3.8 to 12.5 /im point strongly to external heating 
and illumination. Indeed, it seems very likely that 
IRc2 is illuminated and heated by radiation and 
outflow-driven shocks from source I — especially 
if the outflow is oriented along a NW-SE axis 
(Grccnhill ct al. 1998). Assuming that the lumi- 
nosity of source I is a reasonable L ~ 10 4 L Q , 
corresponding to a 10 - 20 M Q star, IRc2 would 
have to intercept ~ 10% of its total energy to 
account for its luminosity. Given the apparent 
relative placement of the IRc2 components and 
source I, this appears to be a plausible number. It 
is possible, however, that both external and inter- 
nal heating play a role in the emission from IRc2, 
as would be expected if a low-mass, dense proto- 
stcllar cocoon were placed very near an energetic 
source. Such objects would be the evolutionary 
forebears of the "proplyds" seen near the Trapez- 
ium. 

It is interesting to note that IRc2 A is clearly 
coincident with the shell H2O masers NE of source 
I, especially at 5.0 /im (Figure 6). If the outflow 
from source I is instead oriented along a NE-SW 
axis, as suggested by Grccnhill et al. (2003), then 
IRc2 A may be thermal emission from the cavity 
being excavated out of the hot core by source I. 
The rest of the IRc2 sources might be shadowed 
by Fs circumstellar disk and hence would be self- 
luminous and/or heated by source n and BN. 
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4.3. The OMC 1 Outflows 

As discussed in Section 1, observational evi- 
dence from a variety of tracers indicates that there 
are (at least) two outflows emerging from OMC 
1. Radio CO and near-IR H2 line emission show 
a fast, poorly collimated bipolar outflow along a 
NW-SE axis. On arcminute scales, the H 2 fin- 
gers appear to be erupting roughly orthogonal to 
the NH 3 ridge (Figure 7). Closer in, the shocked 
regions are correlated, but not coincident, with 
clouds observed at 12 ^tm (Figure 4), suggesting 
that the H 2 emission arises from the shocked sur- 
faces of warm clouds. A plausible explanation is 
that the fingers represent high velocity ejecta mov- 
ing relatively unimpeded into low density regions 
NW and SE of the core, whereas the masers rep- 
resent strong shocks at the locations where the 
ejecta encounter some dense gas in the OMC 1 
ridge. 

Recent VLBA observations of the SiO masers 
have been interpreted as evidence that the disk 
around source I has a minor axis oriented towards 
PA = 45° (Grecnhill ct al. 2003), roughly 90° 
from the previous interpretation by Grecnhill et al. 
(1998). If this is correct, then there is no obvious 
source for the high-velocity "explosive" outflow 
originating from the hot core and aligned NW-SE 
as traced by CO and H 2 emission. But there are 
several alternative possibilities: First, it is possible 
that the disk orientation has flipped by about 90° 
within the last 1,000 years since the production of 
the H2 fingers. Obviously, this hypothesis invokes 
a set of rather peculiar circumstances. Such an 
orientation change may be induced by the close 
passage of another massive star, or the sudden ac- 
cretion of material with radically different angular 
momentum orientation. Second, perhaps one (or 
both) of the sources, I and n, suffered a major 
eruption that was more or less isotropic. As this 
flow ran into the NE-SW oriented ridge of dense 
gas traced by NH 3 observations, it was blocked. 
Some of the ejecta may have been deflected along 
the density gradient to produce a SE-NW flow. 
The problem with this interpretation is that it 
does not explain why the source I disk did not 
cast a mechanical shadow along its plane. Finally, 
it is possible that the recent 7 mm observations of 
source I do not trace dust emission from a disk. 
Rather, this emission traces a dense thermal jet 
oriented SE-NW, exactly along the major axis of 



the high velocity CO outflow and H2 fingers. In 
this latter picture, the distribution of SiO masers 
within an arc second of source I do not trace the 
surface of a disk, but rather the shear layer formed 
where the jet interacts with a thick torus or cored- 
sphere of circumstellar material. While the geom- 
etry of this hypothesis is compelling, the observed 
velocities of the SiO masers and the CO and H 2 
fingers do not match: the SiO masers are geneer- 
ally red-shifted to the NW of source I and blue- 
shifted to the SE while the CO and H 2 emission 
shows exactly the opposite velocity distribution. 

We suggest the following scenario for the 
BN/KL region. Most of the H 2 and OH masers 
(but not the "shell masers" in the immediate vicin- 
ity of source I) are produced by a relatively steady 
bipolar outflow from source n (corresponding to 
the "18 km s" 1 " flow, Genzel et al. 1981) which 
has cleared two cavities to the NE and SW. The 
H 2 maser motions suggest that this outflow has 
been active for at least 3000 years. The relation- 
ship between the explosive outflow associated with 
the H 2 fingers and oriented NW-SE and the out- 
flow from source n is not clear, though both have 
probably contributed to the production of high ve- 
locity gas in the BN/KL region. These flows have 
similar dynamical ages and it is not clear which 
flow stared first. If the explosive ejecta moved 
relatively unimpeded through the medium, then 
the flow from source n may have cleared a cavity 
before the the explosive event traced by the H 2 
fingers. The limited spatial extent of the "shell 
masers" suggests that source I may be younger 
than source n. 

An interesting feature of the NE and SW out- 
flow cavities is the apparent anti-correlation of OH 
and H 2 masers (Figure 5): the NE cavity is dom- 
inated by OH masers and the SW by H 2 masers. 
This discrepancy might be because of the pho- 
todissociation of H 2 in the NE cavity by strong 
UV radiation from source I, n, BN, or perhaps 
even 9 1 C Orionis. The SW cavity may be shielded 
from any UV flux from source I and BN (though 
not n) allowing H 2 to survive and form maser 
spots. 

4.4. IR source "t" 

The mid-IR and H 2 emission north of IR source 
"t" (SE of IRc2; Figure 4) trace a semi-circle of H 2 
emission with a diameter of ~ 3" (1370 AU at 450 
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pc). The 12.5 /jm emission south of IRcll forms a 
foot-shaped ridge following the H 2 emission north 
of source t. In addition, the mid-IR sources IRc8 
and 22 appear to continue the circular structure 
around source t to the south. 

H2 emission can arise in both shock-excited re- 
gions and gas heated to a few thousand K by UV 
radiation. The circular structure around source 
t is highly suggestive of a PDR expanding into 
a higher density cloud - which can be accommo- 
dated if source t is a hot star on the near side, 
and close to the surface, of the hot core. If this 
hypothesis is correct, then source t must also be 
a significant source of UV flux, though not nec- 
essarily in the Lyman continuum since H2 is ex- 
cited by soft-UV photons. Thus, source t may be 
a late B or early A star. As the PDR advances 
outward from source t, the radiation front illumi- 
nates over-dense regions which then emit in the 
mid-IR. We postulate that the local extinction is 
somewhat higher to the south of t, blocking 2 /jm 
H2 emission from the southern part of the PDR, 
while allowing 12.5 /im emission from IRc8 and 22 
to pass. The interface of the PDR with the IRcll 
"foot" seems to be viewed nearly edge-on. 

5. Summary & Conclusions 

In this paper we have presented the highest spa- 
tial resolution map at 12.5 /xm of the BN/KL re- 
gion to date. Despite relatively poor flux calibra- 
tion owing to variable atmospheric transmission 
we were able to detect 6 new sources and resolve 
new morphologies for existing sources leading to 
the following results: 

• IRc2 is resolved into multiple sources, as ob- 
served at ~ 4 ^m(Dougados et al. 1993). 
The point-like nature of sources A - D at 
3.8 fim coupled with the relatively high color 
temperature of IRc2 in the mid-IR, seems to 
suggest that these are self-luminous sources. 
Near-IR polarization maps, however, suggest 
that IRc2 is externally illuminated. Com- 
parison of our 12.5 fim map with previous 
datasets clearly show that the shorter wave- 
length sources all appear at the periphery 
of the 12.5 /xm emission. This lends some 
weight to the picture in which IRc2 is illumi- 
nated and heated by source I, but we cannot 
rule out the possibility that IRc2 is a small 



group of nearby protostars. 

• The elongated structure of near-IR source n 
at 12 fim, coupled with the geometry of the 
mid-IR and NH3 emission, and the location 
of H2O and OH maser spots, is highly sug- 
gestive of a disk with rotation axis along a 
NE-SW axis on the sky. Thus, source n may 
be responsible for the low- velocity 18 km s _1 
flow emerging from OMC 1. 

• In addition, we have discovered an arc of 
emission associated with the BN object and 
some interesting morphology for IR source 
t, SE of IRc2. The BN SW arc is a curious 
feature that requires follow-up imaging and 
spectroscopy at other wavelengths to deter- 
mine its true nature. 

Owing to the extreme obscuration on the line of 
sight, our mid-IR observations are only sampling 
the most massive end of the IMF in BN/KL. The 
stellar density in the region is likely to be very 
high (as suggested by the source I-IRc2 group), 
perhaps even greater than the Trapezium cluster. 
This high stellar density may be due to environ- 
ment or evolutionary state (or both): It is possi- 
ble that the Trapezium Cluster also began with a 
much higher stellar density. Additional deep mid- 
IR observations, especially at 4 and 8 /jm accom- 
panied by 20 /an images to determine spectral en- 
ergy distributions, are required to determine the 
full stellar content of this rich environment. 
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Fig. 1— Filled contour map of the Kcck/LWS 12.5 /^m mosaic for BN/KL. Flux is mJy per 0.081 x 0.081" 
pixel. Contours shown arc for 10, 20, 30, 40, 50, 100, 200, 300, 400, 500, 600, and 6000 mJy; representative 
contours are marked in the figure. The dashed box shows the area not covered by our observations. Coor- 
dinates are centered on the BN object (RA = 05h35ml4.117s, DEC = -05°22'22.90"[J2000]). See Table 1 
for IR source identifications. Continuum radio sources (RS)'T' and "L" are also indicated (Menten & Reid 
1995). 
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Fig. 2. — Filled contour map of the 12.5/itm mosaic in the IRc2 region (grayscale, same contours as Fig. 1) 
overlaid with (a) L' emission and (b) 4.0 (dashed) and 5.0 (solid) /Ltm emission (Dougados et al. 1993) in 
red contours. The spatial resolution of the ~ 4^m data is FWHM = 0.5". See Table 1 for IR source 
identifications. 




Fig. 3. — Detailed 12.5 /xm sub-images (10" x 10") enhancing low level emission around near-IR source n 
(a) and the BN SW arc (b). 
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Fig. 4. — Three-color composite image of BNKL region (31" x 27"): Red = 12.5/im (this work); Green = 
NICMOS 2.15 fim continuum; Blue = 2.12 ^m H 2 emission (NICMOS N212 - N215). A logarithmic stretch 
has been applied to each channel independently. To avoid confusion, stellar features in the H2 channel (blue) 
have been masked out. The NICMOS data were obtained from the Multi-Mission Archive at Space Telescope 
and were first published by Schultz et al. (1999). 
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Fig. 5. — Filled contour map of the 12.5 /zm mosaic (grayscale, same contour levels as Fig 1) overlaid with 
NH3 emission at Vi sr = 6.79 km s _1 (green contours) tracing the Orion "hot core" (Wilson et al. 2000), H2O 
masers (filled dots, Genzel et al. 1981; Gaume et al. 1998), and OH maser clusters (open circles, Johnston 
et al. 1989). Red and blue indicate Doppler shift for the maser(s) relative to Vi sr = 5 km s _1 , the approximate 
velocity centroid of the hot core. The 18 km s _1 H 2 maser outflow center of expansion (COE) is remarkably 
close to IR source "n" . Notice the significant anti-correlation of NH 3 and 12 /im emission, suggesting that 
most of the mid-IR emission comes from at, or behind the hot core. The mid-IR emission and maser spots 
appear to trace bipolar cavities centered on source "n" . Offsets are relative to BN in arcseconds. 
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Fig. 6. — Filled contour map of the 12.5/mi mosaic (grayscale, same contours as Fig. 1) overlaid with NH3 
emission at vi sr — 5.57 km s" 1 (green contours) tracing the Orion "hot core" (Wilson et al. 2000), H 2 
masers (filled dots, Genzel et al. 1981; Gaume et al. 1998). Red and blue indicate Doppler shift for the 
maser(s) relative to vi sr — 5 km s , the approximate velocity centroid of the hot core. The ammonia 
emission strongly suggests a cavity surrounding radio source I with a "thumb" of dense material occulting 
the southern part of IRc2. The "shell masers" surrounding source I may arise in shocked regions in the 
cavity wall. Offsets are relative to BN in arcseconds. 
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Table 1 
Mm-IR Sources in BN/KL 







L^nser 


.Positional 




Aperture 






Name 


RA (") 


Dec. (") 


Error (") 


Flux (Jy) b 


Radius (") 


Notes 


Refs. 


BN 


0.0 


0.0 


0.1 


[600] 


2.025 




2,3 


BN SWarc 


-1.4 


-0.6 


0.1 




Elongated along PA = 131° 


1 


IR n 


3.5 


-9.8 


0.1 


14.2 


0.81 


3,4 


IRc2 A+B 


5.7 


-7.0 


0.1 


[120] 


1.215 


Single 12 (im source between A & B 


5 


IRc2 C 


4.6 


-7.8 


0.1 




6 


IRc2 D 


5.0 


-7.0 


0.1 








6 


IRc2 E 


4.7 


-6.3 


0.1 








1,6 


IRc3N 


-1.6 


-7.2 


0.1 








1 


IRc3S 


-1.4 


-7.9 


0.1 








1 


IRc4 


1.1 


-11.1 


0.1 






Position here is peak at NW vertex. 


3,7 


IRc5 


0.5 


-14.5 


0.2 


19 


0.891 


3,7 


IRc6 


0.4 


-4.5 


0.1 






Central source 


3,8 


IRc6N 


-0.2 


-3.2 


0.1 








3 


IRc6E 


2.0 


-4.8 


0.1 








1 


IRc7 


2.8 


-7.8 


0.1 


108 


1.215 




3,8 


IRc8 


8.4 


-12.2 


0.2 








3,8 


IRcll 


8.8 


-7.9 


0.2 


18.7 


0.81 


Position and flux are for 
primary elliptical source 


3 


IRcl2 


10.3 


-8.4 


0.2 


10.1 


0.81 




3 


IRcl3 


14.5 


-4.8 


0.3 






Diffuse, broken structure 


3 


IRcl4 


14.5 


-7.2 


0.3 






Diffuse, broken structure 


3 


IRcl5 


3.1 


-2.3 


0.1 






Position uncertain (edge) 


3 


IRcl6N 


4.9 


1.3 


0.2 






1 


IRcl6S 


5.5 


-0.8 


0.2 








1 


IRcl6W 


4.0 


0.2 


0.1 








1 


IRcl7 


2.5 


2.8 


0.1 








3 


IRcl8 


6.8 


1.8 


0.2 


2.6 


0.729 




3 


IRcl9 


5.7 


4.2 


0.2 






Associated with NIR "q" ? 


!' 4 


IRc20 


0.7 


-6.9 


0.1 


7.1 


0.729 






IRc21 


4.4 


-5.0 


0.1 


5.2 


0.729 




i 


IRc22 


6.4 


-11.3 


0.2 






Faint N-S ridge 


i 



a Offset from BN (RA = 05h35ml4.117s, DEC = -05°22'22.90", J2000) 

b Total sky-subtracted flux measured in the aperture listed; the sky value was determined from an annulus just beyond the 
aperture. 2a error for total flux is 10% - 27%, see text. Fluxes in brackets (BN and IRc2) are from Gezari et al. (1998) and were 
used for calibration. 

References. — (1) This work; (2) Bccklin & Neugebauer (1967); (3) Gezari et al. (1998); (4) Lonsdale ct al. (1982); (5) Chclli 
ct al. (1984); (6) Dougados ct al. (1993); (7) Riekc et al. (1973); (8) Downes et al. (1981) 



